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Abstract 

We study the singular ordinary differential equation 

dU 1 



dt C{U) 



<t)s{U) + (f>„s{U), (0.1) 



where U € R'^, the functions cps G R'^ and (fins G R^ are of class and is a real valued function. 
The equation is singular in the sense that ({U) can attain the value 0. We focus on the solutions of 
(|0.ip that belong to a small neighbourhood of a point U such that (j>s{U) = (j)ns{U) = 0, C,{U) = 0. We 
investigate the existence of manifolds that are locally invariant for (|0.ip and that contain orbits with a 
suitable prescribed asymptotic behaviour. Under suitable hypotheses on the set {U : C,{U) = 0}, we 
extend to the case of the singular ODE (|0.1|l the definitions of center manifold, center stable manifold 
and of uniformly stable manifold. We prove that the solutions of (lO.lf) lying on each of these manifolds 
are regular: this is not trivial since we provide examples showing that, in general, a solution of (|0.ip is 
not continuously differentiable. Finally, we show a decomposition result for a center stable and for the 
uniformly stable manifold. 

An application of our analysis concerns the study of the viscous profiles with small total variation for a 
class of mixed hyperbolic-parabolic systems in one space variable. Such a class includes the compressible 
Navier Stokes equation. 

Key words: singular ordinary differential equation, stable manifold, center manifold, invariant manifold. 

1 Introduction 

In this work we study the singular ordinary differential equation 

^ </),([/) +0„.(C/). (1.1) 



dt C{U) 



In the previous expression, U G and the functions (j)s and 0„s are (continuously differentiable with 
continuously differentiable derivatives) and take values into M^. The function ( is as well regular and it 
takes real values. We say that the equation is singular because C{U) can attain the value 0. 

Equation (jl.ip is related to a class of problems studied in singular perturbation theory. Consider system 

edx/dt = /(x, y, e) ,^ 2) 

dy/dt = .g(x, y, e), 

where x and y are vector valued functions, £ is a parameter. In singular perturbation theory one is typically 
concerned with the limit e — > and with the corresponding behaviour of the solution (x, y) . Note that (jl.H) 
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can be viewed as as an extension of (|1.2p . in the sense that (|1.2p can be written in the form p.ip : in this 
case, the singularity CiU) in (jl.ip is identically equal to e and hence d(/dt = 0. 

Being the literature concerning (|1.2p extremely wide, it would be difficult to give an overview here. 
Consequently, we just refer to the notes by Jones [TT] and to the rich bibliography contained therein. In 
particular, ^ provides a nice overview of Fenichel's papers [101 EJ [9]. The works (TOllHlHj provide several 
ideas and techniques exploited in the present paper. 

The main novelty of the present work is that we consider the case C is a nontrivial function of the unknown 
U. In particular, this means that dC^/dt ^ in general and hence that we have to face the possibility that 
C(C^(0)) 7^ 0, but Q{U{t)) = for a finite value of t. This is exactly what happens in the examples ()2.14p 
and (|2.19p discussed in Section [2] here. Other examples are provided in a previous work by the same authors 
[5], Section 2. Note that, in all these cases, there is a loss of regularity at the time to at which C((C^(^)) 
reaches the value 0, = minji G [0, +cxd[: C,{U{t)) = O}. More precisely, the first derivative dU/dt either 
has a discontinuity or blows up at t = to. 

Our goal here is to study the solutions of (jl.ip that lie in a neighborhood of a point U such that 
(j)s{U) — 4)ns{U) = and C(^) = 0- We are concerned with the existence of invariant manifolds. More 
precisely, the problem is the following. 

Consider first the case of the non singular ODE 

^ = (1.3) 

and assume that the point U is an equilibrium, namely f{U) = 0. In a neighbourhood of U one can define a 
center and a center stable manifold, which are both locally invariant for (|1.3p . We recall here that, loosely 
speaking, a center stable manifold contains the orbits of (|1.3p that as a; ^ +oo either do not blow up or 
blow up more slowly than e''*, where 77 is a small enough constant depending on the system. More precisely, 
the orbits that lie on a center stable manifold are those having the asymptotic behaviour described before 
and solving a suitable system which, in a small neighbourhood of U , coincides with the original one (11. 3p . 

Also, assume that there exists a manifold E entirely constituted by equilibria of ()1.3p and containing 
U . We are interested in the uniformly stable manifold relative to E. By uniformly stable manifold we 
mean the slaving manifold that contains all the orbits that when t +00 decay with exponential speed to 
some point in E. Note that the uniformly stable manifold does not coincide, in general, with the classical 
stable manifold. Indeed, the stable manifold contains the orbits that decay exponentially fast to the given 
equilibrium U, while on the uniformly stable manifold we only require that the limit belongs to E. 

The existence of a center stable and of the uniformly stable manifold can be obtained as consequence of 
the Hadamard-Perron Theorem, which is discussed for example in the book by Katok and Hasselblatt [13j . 

In the present paper we prove that, under suitable hypotheses, one can extend the definition of center, 
center stable and of uniformly stable manifold from the case of the non singular ODE p.3p to the singular 
case ()l.ip . The manifolds we define are all locally invariant for p.ip and satisfy the following property: 

(P) If U is an orbit lying on the manifold and ({U{0)) ^ 0, then C{U{t)) ^ for every t. 

This, in particular, rules out the losses of regularity (blow up or discontinuity in the first derivative) mentioned 
before. 

We proceed as follows. First, we consider the non singular ODE 

^=0.([/) + C(t/)^„s(C/). (1.4) 

Formally, (|1.4p is obtained from (II. ip via the change of variable r = T{t), defined as the solution of the 
Cauchy problem 

( dr _ 1 

Tt ~ Wm (1.5) 
I ^(0) = 0. 
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However, the function T{t) is well defined only if C[C/(t)] 7^ for every t. In the work we always refer to the 
formulation (jl.4p and we prove the existence of locally invariant manifolds satisfying property (P). We then 
show that a posteriori the change of variable (jl.Sp is well defined and system (|1.4p is equivalent to (jl.ip on 
these manifolds. 
We assume that 

1. For every center manifold of (|1.4p . the intersection between the set {U : C{U) ~ 0} and Ai^ 
contains only equilibria. 

We then define a manifold of slow dynamics as a center manifold of (II. 4p (any center manifold works). To 
simplify the exposition, in the following we fix a manifold of slow dynamics. To define the manifold of fast 
dynamics we assume 

2. there exists a one-dimensional manifold which is transversal to the set {U : ({U) = 0} and is entirely 
constituted by equilibria of (jl.4p . In the following, we denote by E this manifold. 

The manifold of fast dynamics is then defined as the uniformly stable manifold of (jl.ip relative to the 
manifold of equilibria E. Namely, all the fast dynamics converge exponentially fast to some equilibrium in 
E. 

The manifolds of slow and fast dynamics can be regarded as extensions to the general case of the notions 
of slow and fast time scale discussed for example in Fenichel [lOj in relation to system p.2p , namely to the 
case ^ is a parameter. 

We also assume that 

3. The singular set {C(U) — 0} is invariant for both (|1.4p and for the solutions of (jl.ip that lie on the 
manifold of the slow dynamics. 

It is not hard to show that, as a consequence of the previous assumptions, if we restrict system ()1.4p on 
the manifold of slow dynamics, then (jl.4p is equivalent to (|l.ip . We can thus go back to the original variable 
t and get that the solutions of (jl.ip lying on the manifold of the slow dynamics satisfy a non singular ODE. 
We then define a center manifold of (|l.ip as a center manifold of the system reduced on the manifold of slow 
dynamics fTheorem l4.ip . In this way property (P) is automatically satisfied on any center manifold and the 
losses of regularity are ruled out. As before, by loss of regularity we mean the blows up or discontinuities 
in the first derivative that may be exhibited by the solutions of (|l.ip . as shown by the examples (|2.14p and 
(Prg)) in Section H here. 

To extend to the case of the singular ODE (|l.ip the definition of center stable and uniformly stable 
manifold we need some more work. As mentioned before, thanks to assumption 2, there exists a manifold of 
equilibria transversal to the singular surface: we denoted this manifold by E. To define the uniformly stable 
manifold of p.ip relative to E we need to study the solutions of (II. ip which converges to an equilibrium 
on E with exponential speed. Note that this speed can be either bounded or unbounded as C ~* 0; so 
we are looking for a composition of both fast and slow dynamics. Also, to define a center stable manifold 
loosely speaking we have to study orbits that are local solutions of (|l.ip and that do not blow too fast when 
t — > +00. Therefore, we have to deal again with a composition of slow and fast dynamics. 

In both cases (uniformly stable and center stable manifold) the analysis can be seen as an extension of 
the exponential splitting methods for non singular ODEs like ()1.3p . However, as mentioned before what a 
priori can go wrong is that in the change of time scale defined by the Cauchy problem 

' dr _ 1 

, dt'cWW] (1.6) 
^ r(0) = 0. 

some regularity is missing. 

The main result of this paper is the following (a more precise statement is given in Theorem 14. 2p : 
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Theorem 1.1. There is a sufficiently small constant S > such that the following holds. In the ball of 
center U and of radius 6 in M.^ one can define two continuously differ entiahle manifolds and Ai'^^ which 
are both locally invariant for p.ip . The first one, , is the uniformly stable manifold of (|l.ip relative 
to E, while Ai'^^ is a center stable manifold for (11. ip . If U is a solution satisfying C(J7(0)) ^ and lying 
on either A4'^ or M-'^^ , the Cauchy problem (|1.6p defines a diffeomorphism t : [0, +oo[— > [0, +cx)[. In other 
words, if we restrict to either or Ad'^'^ , then the formulations (jl.ip and (|1.4p are equivalent, provided 
that (^([/(O)) 7^ 0. In particular, property (P) is satisfied on both and A4'^^ . 

Also, if U(t) is a solution lying on either Al* or A^'^", then it can be decomposed as 

UiT)^Uf{T) + Usi{r) + UpiT), (1.7) 

where Usi lies on the manifold of slow dynamics and Uf{T) is exponentially decreasing to 0. The perturbation 
term Up{T) is small in the sense that 

|C/^(T)|<fc^|C(C/(0))||C//(0)|e— /4 
for suitable positive constants c, kp > Q. 

From the technical point of view, the key points in the proof of Theorem p.ip are the following two. 
First, we introduce a change of variables which allows us to write system 

^=0.([/) + C(t/)0„.(C/). (1.8) 

in a more convenient form. The precise statement is given in Proposition 14.11 This change of variables 
exploits many of the ideas which in the case of equation (|1.2p lead to the introduction of the so called 
Fenichel Normal Form (see Jones |llj and the references therein) . Here, however, we have to rely on the 
assumptions 1, 2 and 3 discussed above. 

The second main point in the proof of Theorem 11.11 is the analysis of a family of slaving manifolds for 
system (|1.8p . This analysis exploits the presence of a splitting based on exponential decay estimates and 
it is in the spirit of Hadamard Perron Theorem (see for example the book by Katok and Hasselblatt [T5]). 
Here the main results are Theorem 13.11 and Proposition 13.11 Loosely speaking. Proposition 13.11 tells us the 
following. Fix a manifold S, locally invariant for ()1.8p and entirely made by slow dynamics. Then there 
exists a slaving manifold containing orbits that decay to an orbit in S exponentially fast, with respect to 
the T variable. Also, Proposition 13.11 ensures that any solution U lying on the slaving manifold admits a 
decomposition like (|1.7p . namely 

U{T)^UfiT)+U,l{T)+Up{T) 

where Usi lies on S and Uf{T) is exponentially decreasing to 0. The perturbation term Up{T) is small 
and disappears when C{U) = 0, so on the singular surface {U : ((U) = 0} there is no interaction, but 
a complete decoupling. Actually, in the statement of Proposition 13.11 we consider slightly more general 
conditions ensuring that the interaction term disappears. However, the case C{U) = is the one we exploit 
in the following. From the technical point of view, the most complicated point in the analysis is proof of the 
regularity of the slaving manifold, since it involves studying the Frechet differentiability of suitable maps 
between Banach spaces. 



An application of our analysis concerns the study of the viscous profiles with small total variation for a 
class of mixed hyperbolic-parabolic systems in one space dimension. The connection between these viscous 
profiles and the singular ordinary differential equation (jl.ip is discussed in [S] , where we also explain what we 
mean by viscous profiles and by mixed hyperbolic-parabolic systems in this context. In [6j we also discuss a 
remark due to Frederic Rousset [T3] about the Lagrangian and the Eulerian formulation of the Navier Stokes 
equation. Loosely speaking, the connection between viscous profiles and singular ODEs like (|l.ip is that the 
equation satisfied by the viscous profiles may be singular when the system does not satisfy a condition of 
block linear degeneracy defined in f5]. In particular, this happens in the case of the Navier Stokes equation 
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written in Eulcrian coordinates. As we see in Section 12.1.11 the analysis developed in the present paper 
applies to the study of the viscous profiles of the Navier Stokes. 

Viscous profiles provide useful information when studying the parabolic approximation of an hyperbolic 
conservation law. If one restricts to systems with small enough total variation, it is often meaningful to 
focus on viscous profiles lying on a center, a center stable or on the uniformly stable manifold. The literature 
concerning these topics is extremely wide. Here, we just refer to the books by Dafermos |7j and by Serre 
[IB] and to the rich bibliography contained therein for a discussion about the parabolic approximation of 
conservation laws. For the applications of the viscous profiles to the study of this approximation, see for 
example Bianchini and Bressan j5j and Ancona and Bianchini Concerning the analysis of viscous profiles, 
we only refer to Benzoni-Gavage, Rousset, Serre and Zumbrun j^, to Zumbrun [18], and to the references 
therein. For an alternative approach to the analysis of the viscous profiles of the compressible Navier Stokes 
equation, see Wagner [T7] and the references therein. 



The exposition is organized as follows. In Section [Ol we discuss a toy model and we introduce the results 
that are extended in the following sections to the general case. 

In Section [2] we define our hypotheses and in Section r2. 1.1 1 we show that they are satisfied by the viscous 
profiles of the compressible Navier Stokes equation. Also, in Section 12.21 we discuss two examples: each of 
them show that, if one different hypothesis is not satisfied, then the first derivative dU/dt of a solution U of 
(jl.ip may blow up in finite time. 

In Section [3] we define a class of invariant manifolds for an equation with no singularity in it (Theorem [5TT] 
and Proposition [33J- This analysis is applied in Section [H to study the singular ODE (jl.ip . In particular, in 
Section f4. II we define the notions of slow and fast dynamic and we extend the definition of center manifold to 
the case of the singular ODE (jl.ip . In Section 14.21 we discuss how to extend the notions of uniformly stable 
and center stable manifold: the main result here is Theorem l4.2l Finally, Section [131 is devoted to the proof 
of Proposition 14. 11 a technical result which reduces our system to a more convenient form. 



1.1 A toy model 

In this section we discuss a toy model for system (jl.ip . The goal is introducing in a simplified context the 
analysis that is extended in Section [4] to the general case. All the conditions introduced in Section [2] are 
satisfied by the toy model discussed here, except for Hypothesis [2J Hypothesis [5] is a technical condition and 
prescribes that the functions (j)s and (t>ns in (|l.ip satisfy the following: (t)s{U) + ({U)(j)ns{U) is identically 
when U is out of a small enough neighbourhood of the origin. 

Actually, our toy model can be handled with known geometric singular perturbation theory techniques. 
Indeed, in the present section we assume that the function C in p.l|) is just a parameter, namely 



- 0. (1.9) 



dt 

Also, we focus on the case of a linear system: 

^ = ^A,V + Ar.sV VeR^. (1.10) 

In the following we consider only non negative values of t and we focus on the limit ^ — > 0+. The study of 
the limit C — > 0~ does not involve additional difficulties. 
Consider the system 

^^A,v + CAnsV yeM^ (1.11) 

UT 

which is obtained from (|1.10p via the change of variable t — t C,. In the following, we denote by n_ the number 
of eigenvalues of As having strictly negative real part (each of them counted according to its multiplicity) 
and by n+ the number of eigenvalues with strictly positive real part. We denote by no the multiplicity of the 
eigenvalue and, relying on Assumption 1 in the introduction, we assume that there are no purely imaginary 
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eigenvalues. Also, if we write the Jordan form of Ag, then in the block corresponding to the eigenvalue all 
the entries are 0. 

Let C ^ 0+: we are concerned with the behavior of the eigenvalues of the matrix As + C^ns- Thanks to 
results concerning the perturbation of finite-dimensional linear operators (see for example the book by Kato 
|12j . page 64 and foUowings), these eigenvalues can be classified as follows: 

1. n_ eigenvalues converge to the eigenvalues of Ag with strictly negative real part. We denote by M~(C) 
the eigenspace of Ag + (Ans associated to these eigenvalues. 



2. ri+ eigenvalues converge to the eigenvalues of Ag with strictly positive real part. We denote by M+(C) 
the eigenspace of Ag + C^ns associated to these eigenvalues. 



3. the remaining uq eigenvalues converge to as C 
associated to these eigenvalues. 



0+. We denote by M°{() the eigenspace of Ag + (A„ 



When C ^ 0+, the subspace converges to Af^(O), which is the eigenspace of Ag associated to 

eigenvalues with strictly negative real part. The convergence occurs in the following sense: M~(C) is the 
range of a linear application P~(C) G C{M.'^, K.'^). As C ^ 0+, P~{C) converges to P~{0) and the range of 
P~(0) is exactly M~{0). Similarly, when C ^ 0+, the subspaces Af+(C) and Af°(C) converge respectively to 
M^{0) and Af°(0), the eigenspaces of Ag associated to the eigenvalues with strictly positive and zero real 
part. We refer again to Kato [l^ for a complete discussion. 
If V belongs to M~{(), then p. lip is equivalent to 



dV- 



47+CO(l) 



v 



where V~ G M"" and Aj is a n_ x n_-dimensional matrix which does not depend on ( and whose eigenvalues 
have all strictly negative real part. In the previous equation, the entries of the vector V~ are the coordinates 
of u with respect to a basis of A/^(C) and 0(1) denotes a n_ x n_-dimensional matrix which possibly 
depends on ( but remains bounded as C ^ 0+. Its exact expression is not important here. 

If C is sufficiently small, then all the eigenvalues of the matrix [Ag + (0{1)\ have strictly negative real 
part and hence the solution V''{t) converges exponentially fast to 0. More precisely, one has 

|F-(t)| <e--/2|F-(0)|, 

where c > satisfies — c > A for every A eigenvalue of Ag. Coming back to the original variable t, satisfies 

\V-{t)\ < e-=*/2«|l/-(0)|. 



and hence the speed of exponential decay gets faster and faster as C 
as a fast dynamic. 

Conversely, assume that V G Af°(C), then (jl.lip is equivalent to 



O'^. In this sense, we can regard V 



dV^ 
dT 



C 



LoAr^gRo + (0(1] 



where _Ro a-nd Lq are two matrices that do not depend on C. The matrix Rq has dimension N x Uq and 
its columns constitute a basis of Al'^{0). The matrix Lq is no x A'^-dimensional and satisfies LqRq = In^. 
Also, V'^ = LPV and 0(1) denotes an tiq x no-dimensional matrix, which possibly depends on ^ but remains 
bounded as C — *■ 0^. Its exact expression is not relevant here. Coming back to the original variable i, one 
gets 

dV° V ^ 

— ^ LoA„,i?o + CO(l) V\ (1.12) 
dt L J 

and hence V'^ can be regarded as a slow dynamic, because it satisfies the non singular ODE (jl.l2p . 

Applying the same techniques mentioned before, one gets that the eigenvalues of Lo^ns-Ro + can 
be divided into 3 groups: 
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1. eigenvalues that converge to the eigenvalues of LqAjisRo with strictly negative real part. We denote 
by M°~(C) the corresponding cigenspace. 

2. eigenvalues that converge to the eigenvalues of LgAnsRo with strictly positive real part. We denote by 
M°+(C) the corresponding cigenspace. 

3. eigenvalues that converge to the eigenvalues of LoAnsRo with zero part. We denote by M°*'(i^) the 
corresponding cigenspace. 

If V{t) e M°~(C), then V{t) converges exponentially fast to the equilibrium when t +oo, but the speed 
of exponential decay does not blow up as C ^ 0+. 
The space 

Ar(C) = M-iC) ® M°-(C) (1.13) 

can be regarded as uniformly stable space for (jl.lOp . because every orbit entirely contained in this space 
decays exponentially fast to 0. Also, the speed of exponential decay is uniformly bounded from below by a 
constant which does not depend on Another way of interpreting this observation is the following: combine 
equations (|1.9p and (jl.lOp and consider in R'*"'"^ the system 

dC/dt = 

dV/dt = AsV/C + A„,V. ^'-^^^ 

Every point in the subspace {{(, 0)} is then an equilibrium for (|1.14p . Also, for any every orbit V{t) 
lying on M{() converges to the equilibrium (C, 0), and the speed of exponential decay is bounded below by 
a constant independent of 

Conversely, the space M°*'(C) can be regarded as a center manifold of the original equation p.lOp . 

In Section |4] we will extend the previous considerations to the case of the non linear equation 

:|._L_,.(f/)+,,„(c/), 

where C(f^) is in general a non constant function. 



2 Hypotheses 

In this section we discuss the hypotheses assumed in the work. 

More precisely, in Section [TT] we state the conditions imposed on the singular ODE 

f-^)MU) + <PUu). (2.1) 

These conditions can be divided into two groups: Hypotheses [U [2 [21 allow to avoid some technical compli- 
cations, but could be actually omitted at the price of much heavier notations. On the other side. Hypotheses 
El El m [71 [8] are much more important and they will be deeply exploited in Section [H Note, however, that in 
Section [31 we are not directly concerned with the singular ODE (|2.ip and that we do not exploit Hypotheses 

Moreover, in Section we discuss three counterexamples. They show that, if either Hypothesis [7] or 
Hypothesis [HI is violated, then the results discussed in the following sections do not hold. In particular, there 
might be solutions of (j2.ip that are not continuously differentiable. 

Finally, in Section 1.1 1 we verify that the conditions introduced in Section [TT] are satisfied by the viscous 
profiles of the compressible Navier Stokes equation written in Eulcrian coordinates. 



7 



2.1 Hypotheses satisfied by the singular O.D.E 

Define 

F{U) = MU) + C{U)^ns{U) (2.2) 
and consider the non singular ordinary differential equation 

^^F{U) C/eM^. (2.3) 
dr 

Formally, (|2.3p is obtained from (|2.f p via the change of variables r — T{t) defined as the solution of the 
Cauchy problem 

( dr _ 1 

M ~ Q[U{t)] ^2.4) 
I ^(0) = 0. 

However, the function r(i) is well defined only if C[J7(t)] ^ for every t. To overcome this difficulty we 
then proceed as follows: we state all the hypotheses referring to the formulation ()2.3p . Relying on these 
hypotheses, in Section |4] we prove the existence of various locally invariant manifolds for (|2.3p satisfying the 
following property. If U is an orbit lying on one of these manifolds and C[t^(0)] ^ 0, then C,[U{t)] ^ for 
every t. If we restrict to the orbits lying on these manifolds, (|2.3p is equivalent to (|2.ip . 
To simplify the exposition, we assume the following: 

Hypothesis 1. The initial datum [/(O) of satisfies C(t^(0)) > 0. 

The case C(C^(0)) < does not involve additional difficulties. The main difference is that, if C(C/(0)) < 0, 
then the change of variable defined by (|2.4p has negative derivative. As a consequence, when t +oo the 
function r(t) —oo. Loosely speaking, the statements given in the present paper can be extended to the 
case C(C/(0)) < in the following way. All the statements concerning the fast dynamics and referring to the 
stable space or to stable-like manifolds have to be replaced by analogous statements concerning the unstable 
space or unstable-like manifolds. However, we will not consider the case C(f^(0)) < exphcitly. 

Before stating the other hypotheses, we recall that we want to study (|2.ip and (|2.3p in the neighbourhood 
of an equilibrium point U such that F(U) = and C{U) = 0. It is not restrictive to take U — 0. Namely, in 
the following we assume 

F(0) = C(0) = 0. (2.5) 

Also, we can assume the following. Fix a positive constant 5 > and consider a smooth cut-off function 
p{U) satisfying 

( u \u\<s 

P^^^ = I \U\> 25. 
In the following, instead of studying system (|2.3p we focus on 



'^^p{U)F{U). 

However, to simplify the notations instead of writing each time p{U)F{U) we assume that Hypothesis [2] 
holds. 

Hypothesis 2. The function F satisfies the following condition: if \U\ > 26 then F{U) — 0. 
The exact size of the constant 6 will be discussed in the following. 

Note that Hypothesis [2] is not restrictive if the goal is to study the solutions of (|2.3p that remain confined 
in a neighbourhood of the origin of size S. Loosely speaking, the analysis developed in Sections [3] and U] can be 
extended to the orbits of systems that violate Hypothesis [H as far as these orbits remain in a neighbourhood 
of the origin with size d. In particular, the manifold described in Sections [3] and [3] are no more invariant if 
Hypothesis [2] is violated: they are just locally invariant. 
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We also introduce the following simplification. It is not restrictive to assume that all the eigenvalues of 
DF(p) have non positive real part. Indeed, this condition is satisfied if we restrict to a center-stable manifold 
for (|2.3p . As mentioned in the introduction, the existence of a center stable manifold can be obtained as a 
consequence of the Hadamard Perron theorem, which is discussed for example in the book by Katok and 
Hasselblatt [13j (Chapter 6, page 242). Also, note that if ((U{0)) < then it is not restrictive to assume 
that all the eigenvalues of DF{0) have non negative real part: this can be obtained considering the solutions 
that lie on a center unstable manifold. 

Hypothesis 3. The Jacobian DF{0) admits only eigenvalues with non positive real part. 

Also, we assume the following non degeneracy condition: 
Hypothesis 4. The gradient VC(0) ^ 0. 

Let S be the singular set 

S:^{U: C{U)^0}. (2.6) 

Thanks to the implicit function theorem. Hypothesis |4] ensures that in a small enough neighbourhood of 
the set S is actually a manifold of dimension — 1, where N is the dimension of U. 



Hypothesis 5. Let A^^ be any center manifold for (|2.3p around the equilibrium point 0. If \U\ < 6 and U 
belongs to the intersection Ai'^ fl S , then U is an equilibrium for (|2.3p . namely F{U) = . 

Concerning equilibria, we also assume the following 

Hypothesis 6. There exists a manifold of equilibria Ai'^'^ for (|2.3p which contains and which is transversal 
to S. 

Let Ueq be the dimension of Ai^'^. We recall that the manifolds S and Ai'^'^ are transversal if the 
intersection S D Ai'^'^ is a manifold with dimension n^q — 1 (as pointed out before, the dimension of S is 
N -I). 

Hypothesis 7. For every U & S, 

VC{U) ■ F{U) = 0. (2.7) 
Thanks to Hypothesis [7] and to the regularity of the functions C and F, the function 

_ VC(^) • F{U) 

can be extended and defined by continuity on the surface S. 

Hypothesis 8. Let U E S be an equilibrium for (^3]) . namely C{U) = and F{U) = 0. Then 

G{U) = 0. (2.8) 

In Section [TTT] we introduced, in the case of a toy model, the notion of slow and fast dynamics. These 
notions will be extended in Section BTT] to the general non linear case. Hypotheses [5] and [7] can be then refor- 
mulated saying that the set S is invariant for the manifold of the slow and of the fast dynamics respectively. 

Remark 2.1. Consider system 

^ = ^0,({/) +</.„,([/). (2.9) 

Also, let f{U) be a regular, real valued function such that /(O) > 0. Clearly, ()2.9|) is equivalent to 

dU 1 



dt aU)f{U) 



MU)f{U)+(bns{U) (2.10) 



and C{U)f{U) 0+ if and only if C{U) — > 0+, at least in a sufficiently small neighbourhood of U — 0. By 
direct check, one can verify that Hypotheses [T] . . . |H]are verified by the couple {(, F) if and only if they are 
verified by the couple (C/, Ff). 
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Remark 2.2. As we will see in Section 21 Hypothesis [S] can be reformulated saying that the slow dynamics 
intersecting the singular manifold {U : C{U) = 0} are equilibria for system (|2.3p . Heuristically, this means 
that we require that the limit as C(?7(0)) ^ 0+ of a solution of (|2.ip is a solution of the limit system. In 
other words, we want to rule out the possibility of a relaxation effect. 

As it shown by the examples discussed in next section, the assumptions on the invariance of the manifold 
{U : C{U) — 0} with respect to the slow and the fast dynamics (Hypothesis [7] and [5] respectively) are due 
to the fact that we want to have a smooth invertible time rescaling t — t{T) defined by (|2.4p . 



2.1.1 The case of the compressible Navier Stokes in Eulerian coordinates 

In this section we show that Hypotheses G] El ... [HI are satisfied by the ODE for the viscous profiles of the 
compressible Navier Stokes equation written in Eulerian coordinates. Also, Hypothesis |2l is not restrictive if 
the goal is to study the viscous profiles entirely contained in a small neighbourhood of an equilibrium point. 

The case of the Navier Stokes written in Lagrangian coordinates was already discussed for example in 
Rousset [in]. When the equation is formulated using Lagrangian coordinates, the ODE satisfied by the 
viscous profiles is not singular. 

The compressible Navier Stokes written in Eulerian coordinates is 



Pt 







ri 



+ pe+ p 



(2.11) 



VVVr, 



Here, the unknowns are p{t, x), v{t, x) and 9{t, x). The function p represents the density of the fiuid, v is 
the velocity of the particles in the fiuid and 9 is the absolute temperature. The function p = p{p, 0) > is 
the pressure and satisfies Pp > 0, while e represent the internal energy. In the case of a polytropic gas, the 
following relation holds: e — OR/{'y — 1), i? being the universal gas constant and 7 > 1 a constant specific 
of the gas. Finally, i^{p) > and k{p) > represent the viscosity and the heat conduction coefficients 
respectively. 

After some manipulations (see [6J for details), one gets that the equation satisfied by the steady solutions 
of the compressible Navier Stokes can be written in the form 



provided that U — ^p, v, 9, , C,{U) — v and 

F{U) = 



Al^z/aii 

V z 

A22V - A2iA2i/aii 



J 



(2.12) 



The equation satisfied by the travelling waves of the compressible Navier Stokes equation in one space 
variable is similar, the only difference being that the singular value is f = cr, where a is the speed of the 
travelling wave. 

In (|2.12p . A21 is a vector in and A21 denotes its transpose. Also, w = px ad z = (vx, dx) ■ The 
function an is real valued and strictly positive if p is bounded away from 0. The matrix b has dimension 
2x2 and all its eigenvalues have strictly positive real part. The exact expression of these terms is not 
important here. Finally, 

1 / pv - V px pe 

Pe~vvxl9 pvee/O ~ k' px/9 



A22 = 



(2.13) 



Note that A22 depends on px but, plugging w = —A^^z/{aiiv) into (|2.13p one gets that A22V evaluated at 
a point (p, u = 0, 0, z = 0) is the null matrix. 
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Thus, the Jacobian DF satisfies 



/ 



DF{p, v = 0,9,z = 6) 



V -b-^A2iA^Jai, J 



where O2 denotes the 2x2 null matrix. Since ^2i^2i/'^ii admits only eigenvalues with strictly positive real 
part, then DF admits only eigenvalues with non positive real part and hence Hypothesis |3] is satisfied. Also, 
the dimension of every center manifold of 

is 3. Since the subspace {z = 0} is entirely constituted by equilibria of the equation, it coincides with 
the center manifold. Thus, Hypothesis \E\ is satisfied. Since ({U) = v, then Hypothesis [H is also verified. 
Concerning Hypothesis [6l this is satisfied because {z = 0} is transversal to the singular surface {v = 0}. 
Finally, by direct check one can show that also Hypotheses [7] and [5] are verified. Thus the machinery 
developed in this work applies to the study of viscous profiles with small total variation of the compressible 
Navier Stokes equation written in Eulerian coordinates. 



2.2 Examples 
2.2.1 Example 

Example (|2.14p deals with a system which satisfies Hypotheses [U [3] . . . [6] and Hypothesis [H but does not 
satisfy Hypothesis [71 We exhibit a solution of this system which has a blow up in the first derivative and 
hence it is not continuously differentiable. The loss of regularity experienced in Example (|2.14p regards a 
solution U such that C[t^(0)] 7^ 0, but C{U) reaches the value for a finite value of t. 
Consider the following system: 

dui/dt — —U2/U1 

dU2/dt ~ — 7i2 

T 



(2.14) 



Let U = {ui, , C,{U) = ui and 



System (j2.14p can then be written in the form 

dU 1 



dt C{U) 

In this case, the function F{U) defined by ()2.2p is 

FiU) = 



-U2 
-U2U1 



By direct check, one can verify that Hypotheses [T] . . . [5] and Hypothesis [S] are satisfied by (|2.14p . On the 
other side. Hypothesis [7] is not verified in this case. Indeed, the singular surface S defined by (|2.6p is in this 
case the line {ui = 0} and 

VC-F = -U2 

is in general different from on 5. 
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The solution of (|2.14p can be explicitly computed and it is given by 

ui{t) = ./ui(0) + U2(0)(e-*-l) 

(2.15) 

U2(<) =M2(0)e-* 

Choosing ^2(0) > ui(0) > 0, one has that the solution ui{t) can reach the value for a finite t. Note that 
at that point t the first derivative dui/dt blows up: thus, the solution (|2.15p of (|2.14p is not . 

2.2.2 Example ((TTCl) 

Example ()2.2.2p deals with system (|2.16p . which is apparently very similar to ()2.14p . However, in the case of 
([2?T6l) Hypotheses [HIS]. . . Hare all verified. We show the solutions of ((27T61) are regular. Also, if C[C^(0)] 7^ 
then C,[U{t)] ^ for all values of t. 
Consider system 

■ /^"iM = --; (2.16) 

au2/d,t ~ —U2/U1 

Set U = (^ui, U2^ , C(^) = "1- System (|2.16p can be written in the form 

^ 0,(C/) + 0„.((7) 



dt C{U) 
provided that 



MU)-[ ) <t>ns{U)^i I' 



Also, the function F{U) defined by (|2.2p is in this case 

F{U)- 



-U2U1 

-U2 

By direct check, one can verify that Hypotheses [T] . . . [8] are all verified in this case. 
To study system (|2.16p we can proceed as follows. From (|2.16l) we have 

dui/dt U2 , 
= ~ du2 / dt 

Ui Ui 

and hence 



In 

Eventually, we obtain 



ui{t) 
wi(0) 



= U2(0 -W2(0). 



ui(t) =Mi(0)e"^W-"^(°). (2.17) 

Choose Mi(0) > 0. To prove that ui{t) ^ for all t it is enough to show that U2{t) is well defined (and in 
particular finite) for every t > 0. In the following we also prove that U2{t) is also for every t > 0. This 
guarantees that no loss of regularity occurs. 

Plugging (|2.17p into the second line of (|2.16p we get 

du2ldt ^e"=(°)-"^(*). (2.18) 

ui(0) 

Note that W2 = is an equilibrium for (|2.18p . Also, if ^2(0) < then du2/dt > and hence U2{0) < U2{t) < 
for every t. Conversely, if ^2(0) > then du2/dt < and hence < U2{t) < ^2(0) for every t. In both cases, 
we get that U2{t) is well defined and regular for every t > 0. 
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2.2.3 Example ([^1^ 



With Example (|2.2.3p we discuss a system which satisfies Hypotheses[Tl|3]. . . but does not satisfy Hypoth- 
esis [HI As in Example ()2.14p , we exhibit a solution of this system which is not continuously differentiable 
and the loss of regularity regards a solution U such that C[f^(0)] 7^ 0, but ({U) reaches the value for a 
finite value of t. 



Consider the following system: 



dui/dt = — U3 

du2/dt = -U2/U1 (2.19) 
du^/dt ~ — U3 



Let U = 



{ui, U2, , C{U) = ui and 




MU) = I -U2 I (l^nsiU) - I 

-U3 



System (|2l4l) can then be written in the form 

dU _ _J_ 

and the function F{U) defined by (|2.2p is 



-U3U1 

F{U) - I -U2 

-U3UI 

By direct check, one can verify that Hypotheses [T] . . . [7] are verified by ()2.19p . On the other side, Hypothesis 
HI is not satisfied in this case. Indeed, the surface S = {U : C{U) = 0} is the plane {ui = 0}. Thus, the set 
of points such that C{U) = and F{U) = is {ui = U2 ~ 0} and 

^C-DF-(vcY ^-U3 

is in general different from zero on this line. 

An explicit solution of (|2.19p can be obtained as follows. From the third and the first equation we get 
respectively 

usit) = M3(0)e-* 

=Ml(0)-M3(0)+U3(0)e-*. 



Assume that 1/3(0) = Aui{Q) for some constant A whose exact value is determined in the following. The 
equation satisfied by U2 becomes 

dU2 U2 

~dt 



Thus, we obtain 



and hence 



d_ 

di 



In 



Aui(0)e-* + ui(0)(l-A)' 

In (ui(0)(l-A)e* + Aui(0) 



1 



ui(0)(yl- 1) dt 



U2(t) {l-A)e*+A 



(A-l)Mi(O) 



for a suitable constant B. If {A — 1)mi(0) > 1, then the first derivative du2/dt blows up at t = ln(yl/A — 1). 
Note that this is exactly the value of t at which ui{t) attains 0. 

In general, for every wi(0) > if 1/(A — l)wi(O) is not a natural number, then the solution is not in 
C™ for m — [1/{A — 1)mi(0)] + 1. Here [1/{A — 1)^1 (0)] denotes the entire part. Thus, we have a loss of 
regularity in higher derivatives. 
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3 Uniformly stable manifolds 



In this section we consider the system 

^=F([/), (3.1) 

where U £ and the F : is a regular function. We are interested in the behavior of the solutions 

in a small enough neighbourhood of an equilibrium point. We can then assume that such an equilibrium 
point is 0, namely F{0) = 0. Also, we assume 

\U\ > 2(5 =^ F{U) = 0. (3.2) 

Because of Hypothesis [51 (|3.2I) is not restrictive in view of the applications discussed in Section U] Also, 
because of Hypothesis [3] we assume that all the eigenvalues of the Jacobian DF{0) have non positive real 
part. Note, however, that in this Section we exploit none among Hypotheses [H [H O El [3 and [51 



3.1 Notations and preliminary results 

3.1.1 Frechet differentiability of the fixed point of a family of maps 

In the following we have to study the regularity of the fixed points of a family of maps depending on a 
parameter. To do this, we exploit Lemma |3. II Note that the hypotheses there are not sharp and the result 
could be improved. However, to avoid technical complications we restrict to those hypotheses since they are 
satisfied in the cases we discuss in the following. 

Let X be a closed subset with non empty interior in a Banach space X and let Y be an open subset of 
another Banach space Y. Also, let 

T -.X xY X 

be a map such that, for every y G Y, T(-, y) takes values in X and is a strict contraction, namely there 
exists some constant fc < 1 such that 

\\T{xi, y) - T{x2, y)\\x < k\\xi - X2\\x Vxi, X2 G X. 

Thanks to the Contraction Mapping Theorem, we can define a map 

Y ^X 



which associates to y the fixed point of the map T(-, y). We denote this map by x{y) and we are interested 
in its regularity. Assume that, for every y, x{y) belongs to the inner part of X. Also, assume that, for 
every (x, y) d X x Y such that x is an inner point, T{-, y) is Frechet differentiable at x and denote by 
Tx{x, y) G C{X, X) its differential. Also, assume that, for every {x, y) in the same conditions as before, the 
map T{x, •) is Frechet differentiable at y and denote by Ty{x, y) 6 CiY , X) its differential. 

The proof of the following result relies on standard techniques (see for example the book by Ambrosetti 
and Prodi [IJ) and will be therefore omitted. 

Lemma 3.1. Assume that the map x{y) is Lipschitz continuous and fix a point (x, y) such that x = x{y). 
Also, assume that T is Frechet differentiable with respect to the variable y at the point (x, y), namely Ty{x, y) 
exists. IfTx{x, y) is defined and continuous in a neighbourhood of {x, y), then x{y) is Frechet differentiable 
at y and the differential is 



I -Tx{x{y),y) oTy{x{y),y), 



(3.3) 



where I denotes the identity. 



Note that the map 



I — Tx \x{y), y] is invertible because r(-, y) is a strict contraction on X. 
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Remark 3.1. We want to give a sufficient condition to have that x{y) is Lipschitz continous. Assume that 
there exists a constant L such that, for every j/i, y2 ^ 

\\T{x{y,), yi)~T{x{yi), 2/2) |U < i||2/i - ^2 ||y • (3.4) 
Then the map x{y) is Lipschitz continuous. Indeed, 

\\x{yi) - x{y2)\\x = \\T{x{yi), yi) - T{x{y2), 2/2) ||x 

< \\T{x{yi), yi) - T{x{yi), 2/2) \\x + \\T{x{yi), 2/2) - T{x{y2), ^2) ||x 

< L\\yi - y2\\Y + k\\x{yi) - x{y2)\\x 

Since fc < 1, we get that x{y) is Lipschitz continuous. 



3.1.2 First change of variables 

Consider system (|3.ip . Let V~ be the eigenspace of the Jacobian DF{0) associated to eigenvalues with 
strictly negative real part. Also, let V'^ be the eigenspace associated to the eigenvalues with real part. 
Also, fix V", a center manifold of (|3.ip around the equilibrium 0. Finally, let V~ be the stable manifold. The 
manifolds V" and V~ are tangent at the origin to V'^ and V~ respectively. Note that M.^ = V'^ ® because 
DF{0) admits only eigenvalues with non positive real part. Thanks to the local invertibility theorem, in a 
sufhciently small neighbourhood of the origin we can define a local diffeomorphism T such that the following 
conditions are satisfied. Let U = T{U), then U satisfies 

^-nn (3.5) 

where f{U) = DT{T-^{U))F{T-^{U)) . Write U = (X- , X°), where X° has the same dimension as V° 
and X^ has the same dimension as V~ . Then the stable manifold of (|3.5p is the subspace {X" = 0}, while 
the center manifold is the subspace {X~^ = 0}. In the following, we assume that the constant 5 in (j3.ip is 
small enough to have that the local diffeomorphism T is defined in the ball of radius 26 and center at the 
origin. Also, to simplify the notations we do not write U, X~ and but just U, X'^ and X". 



3.1.3 A priori estimates 

We rewrite system p.Sp as 



dx-/dT = f-{x-.x^) 
dxydT^ fix-, x°) 

The subspaces {X~ — 0} and {X° = 0} are locally invariant for (|3.6p since they represent respectively a 
center and the stable manifold. Thus, /~(0, X") = for every X'^ and f{X~, 0) = for every X~ . As a 
consequence, 

f-{X-, X") ^ A-{X-, X°)X- fix-, X") = A°{X-, X^)X° (3.7) 

for a suitable matrices A~ and A^. By construction, A~(0, 0) admits only eigenvalues with strictly negative 
real part and A°(0, 0) has only eigenvalues with zero real part. As a consequence, the following holds. Let 
n_ denote the dimension of X~ and fix a constant c > satisfying ReX < — c for every A eigenvalue of 
A~(0, 0). Then there exists a constant C_ > such that 

VX- eR"-, |e^"(°'°)*X-| < C-e-'*\X-\- (3.8) 
Also, if S is small enough and |X^(0)| < 5, then the solution of the Cauchy problem 

dx-/dT = f-ix-, 0) 

X-(r = 0) = X-iO) 
satisfies 

|X-(t)| <C_e-=^/2|X-(0)|, 
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where c > is as before a constant such that ReX < — c for every A eigenvalue of A (0, 0). 
Plugging (|3.7p in (|3.6p we get 

dX-/dT = A-{X-,X°)X- 

dX^dr = XO)XO. ^ 

In view of the apphcations discussed in Section [3] it is convenient to take into account the following 
situation. Assume that there exists a continuously differentiable manifold Zq containing the stable manifold 
= 0} and satisfying 

f{X-,X°)=0 y {X- , X°) Zq. (3.10) 

Actually, this assumption is not restrictive, in the sense explained in Remark l3.2l at the end of Section [3. 1.31 
Applying, if needed, a local diffeomorphism, we can assume that X° = (C, uq) and that Zq ^ {C ^ 0}. 

Since the stable manifold is entirely contained in Zq, such a diffeomorphism does not produce any change 

on X~ , but only on X'~' . In the following we assume that the constant 6 in Hypothesis O is so small that the 

local diffeomorphism is defined in the ball of radius 26 and center at the origin. 

Consider the system restricted on the center manifold {X^ = 0}: since the subspace {C = 0} is entirely 

made by equilibria, then we get that the equation 

dX°/dT = i°(0, X°)X° 

becomes 

" dC/dT = B(O^C, "o)C (3^^^) 
duo/dr = (7(0, C, "0)^ 

where B and C are suitable matrices. Note that, by construction, B{0, 0, 0) admits only eigenvalues with 
zero real part. Fix a constant e such that ReX < — e < for any A eigenvalue of A~{0, 0): also, we impose 
e < c, where c is the same as in (|3.8p . Assuming that the constant S in Hypothesis [2] is sufhciently small we 
can assume that every solution ( of p. lip satisfies 

|C(r)|<O(l)e-H|^(0)| (3.12) 

for some suitable constant 0(1). Since in (|3.1ip the matrix C is uniformly bounded, we get that 

|"o(t)-wo(0)| <O(l)e^l"l|C(0)| (3.13) 

for a constant 0(1) (possibly different from the one in p.l2p ). We introduce the following notation: given 
a point X^ = {(, uq) on the center manifold we call the point 

y° = (0, uo). (3.14) 

Clearly depends on X^, but to simplify the notations we won't express this dependence explicitly. 
Combining (|3.12p and (I3.13P we then obtain 

|XO(r)-r°(0)| <fcoe^l^l|C(0)| (3.15) 

for a suitable constant fco- 

Finally, note that, since both A^ and A are zero when \{X^, X'^)\ > 26, then any non constant solution 
of (|3.9p satisfies 

\X°{t)\<2S \X-{t)\<26 Vt. (3.16) 

Remark 3.2. The hypothesis that there exists a manifold of zeroes Zq is not restrictive. Indeed, assume 
that the set of the zeroes of coincides with the stable manifold {X" ~ 0}. In this case, we can set ( ~ 
the component uo disappears and given X^ the element is just itself. This notation ensures that the 
estimate (I3.15P still holds. As it will be clear from the the following, the only fact about Zq we exploit in 
the proof of Proposition 13.11 is estimate p.l5p . As a consequence. Proposition 13.11 can be extended to the 
case Zq is just the stable manifold. 

In other words, the presence of a manifold of zeroes wider then the stable manifold is not strictly necessary 
for the existence of the uniformly stable manifold introduced in Theorem 13.11 However, it allows to get a 
sharper estimate in (|3.24p . 
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3.1.4 Linear change of variables 

In the statement of the foUowing lemma we denote by Uc the dimension of X'^, then N — Uc + n_. The 
proof is standard, so we omit it. 

Lemma 3.2. For every M > 0, there exists a linear change of variables W^" —>■ M."'" such that in the new 
coordinates satisfies 

dX^/dr = A{X-, X")X°, (3.17) 

for a suitable matrix such that 

A\6,6)^A° + N°, (3.18) 
where A'-' and N'^ enjoy the following properties: 

|e'^°*X"| < \X°\ yt >0, X" e (3.19) 

and 

\N"X"\ < X° e M"^ (3.20) 

We specify in the foUowing how we chose the constant M. 

Remark 3.3. If we apply the linear change of coordinates introduced in Lemma [3^ then it is no more true 
that X° = {(, uq) where {( = 0} is the manifold Zq of equilbria for However, estimate p.lSp still holds, 
provided that we change if needed the value of the constant fcp and we take, instead of X°(r), Y_q and C(0), 
their images trough the linear change of variables. 

3.2 Uniformly stable manifold of an orbit 

We are now ready to introduce Theorem l3.1l In formula (|3.24p . C is the component of X'^ ~ {(,, uq) according 
to the decomposition introduced in Section [3. 1.21 

Theorem 3.1. Let Hypotheses\M and\3[hold. If the constant 5 in Hypothesis\^is sufficiently small, then the 
following holds. 

Fix an orbit Y°{t) = (O, X^{t)) of 

dX-/dT^A-{X-,X")X- 

that lies on the center manifold and satisfies \X^(0)\ < S. Then we can define a uniformly stable manifold 
relative to Y'^ (r) . This manifold is defined in the ball of radius 6 and center at the origin, is parameterized by 
{X'^ = 0} and is tangent to this subspace at the origin. Also, it is locally invariant for p.21|) . meaning that 
if the initial datum lies on the manifold, then (X~(r), X'^(t)) belongs to the uniformly stable manifold for 
\t\ sufficiently small. Every orbit lying on the uniformly stable manifold relative to Y'^{t) can be decomposed 
as 

X{t)^Y"{t)+Y-{t) + UP{t), (3.22) 
where the components Y^ = (X^(t),0) and UP{t) satisfy respectively 

|X-(r)| < fc_|X-(0)|e-'=^/2 (3.23) 

and 

|C/p(t)| <fcp|C(0)||X-(0)|e-^-/4. (3.24) 

In ()3.23p and (|3.24|) . c, fc_ and kp are suitable constants. In particular, c is the same as in p.Sp . 
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3.3 Proof of Theorem 

Let the orbit (0, X^{t)) be given. We denote by — X^{0) and by the corresponding projection, 
defined as in (|3.14p . Also, if we write — (C(0), uo(0)) then we set 

C = C(0). (3.25) 

By definition, X'^{t) is a solution of the Cauchy problem 

dx7dT = io(o, x°) 



The proof of Theorem 13. II is divided in several steps: in Section 13.3.11 we introduce the spaces of functions 
we exploit in the proof. In Section r3. 3. 21 we are concerned with the component Y^{t) in p.22p . In Section 
13.3.31 we deal with the "perturbation" term UP{t) in (|3.22p . Both the components Y~{t) and U^ir) are 
obtained as fixed points of suitable contractions: in Section [3.3.4l we study their regularity applying Lemma 
13.11 Finally, in Section 13.3.51 we conclude the proof of Theorem 13.11 putting together all the considerations 
carried on in Sections 13.3.11 13.3.21 13.3.31 and 13.3.41 



3.3.1 Definition of the functional spaces 

Let n_ denote the dimension of X^ . Also, ric denotes the dimension of as in the statement of Lemma 
In the following we exploit the following Banach spaces of functions: 

y- = [X- e C"([0, +oo[, R"- ) : < +00} (3.27) 

and 

y = [x° e C°([0, +oo[, ) : \\X°\\o < +00}, (3.28) 
where the norms || • ||_ and || • ||o are defined as follows: 

||X-||_=sup{e-/2|X-(r)|} ||XO||o = sup{e-^H|x"(r)|} (3.29) 

r r 

The constants c and e are the same as in p.8p and (|3.15p respectively. Also, we consider the following closed 
subsets of and Y^: 

ys = {X- e C"([0, +oo[, R"- ) : < k^S.} y^ = {x" e C°([0, +oo[, R"^ ) : \\XX < ko6.] 

(3.30) 

We specify in the following how to determine the exact value of the constant fc_, while the constant fco is 
the same as in p.l5|) . Also, note that the spaces and y^ are equipped with the same norms || • ||_ and 
II • II as y" and y° respectively. 

We will also need the space of functions defined as follows. Let c be as p.8p and let a £ [0, c[. Consider 
the space 

yP = {([/-, U°) e C°([0, +oo[, R"=+"- ) : \\X-\\p,rt < +00} (3.31) 
which depends on a because it is equipped with the norm 

Y 

Also, we will exploit the closed subset 

y'L - {{U-, t/") e C"([0, +oo[, R"=+"^ ) : \\{U-, f/°)||pert < kpS'}, (3.32) 

which is equipped with the same norm as y^. We specify in the following how to determine the values of 
the constants k„ and a. 



\\{U-, U°)\\pert=snp\e 



{c+a)T/4 



|(C/-(r)| + |C/°(r)|) 
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3.3.2 Analysis of the stable component 

This step is devoted to the definition of Y^{t) = {X^{t), 0). Fix a vector G IR"" satisfying < S. 
We define X~(r) as the solution of the Cauchy problem 



dX-/dT = A-{X-, Y")X- 

x-{o) = x-, 



(3.33) 



where is given by p.l4p . It is known that, for any fixed and X , X can be obtained as the fixed 
point of the application 

T- : yi ^ 

defined by 



T-{X-)[t] = e'^'^X- 



A-{X-{s), r°) - A- X-{s)ds 



(3.34) 



where A~ = ^"(0, 0). The space is defined in p.30p . More precisely, if the constant /c_ in p.30p satisfies 
k- < C- and the constant 6 is p.30p is sufficiently small, then the map T~ takes values in y^ and is indeed 
a contraction. Also, the fixed point satisfies 



(3.35) 



We are now interested in the differentiability of the fixed point with respect to Y^ and A . To study it, we 
recall that 

2o = {(Z",0,iio)}CM^ 

We then regard T~ as an application 

:Zox y- -> y- 

and we verify that the hypotheses of Lemma [3.11 are satisfied. The space y^ is defined by p.27p . The 
Frechet derivative of with respect to (A~, Y^) is a linear map e ^{Zq, y^)- Evaluated at the point 
{h~ , li) & Zq it takes the value 



r-(/i-,/i")[T] 



T J 



L>yoA"(A"(s), i:°)[/i°] X-{s)ds 



In the previous expression, DyoA (A (s), 11°) [/i''] denotes the differential of the function A (A (s), 
with respect Y^, applied to the vector h^. If A~ = 6 then, no matter what Y_^^ is, the differential maps 
{h~ , ll') into the function "^bT . 

The Frechet derivative of T~ with respect to A^ is a linear map G C{y^ , y^)- Evaluated at the 
point G it takes the value 



S-{h-)[T] = 



^^-^^{[A-{X-is),Y°)^A-]h-{s)+[D^-A-{X-{s),Y°)[h-{s)]\x-{s)} 



>ds 



In the previous expression, Dx~A-{X-{s), Y") [h-{s)] denotes the differential of the function A" (A-(s), Y°) 
with respect A~, applied to the vector h~[s). 

Both T~ and S~ are continuous if viewed as maps from Zq xj^a to C{Zo, y~) and C{y^ , y~) respectively. 
Thus, the hypotheses of Lemma |3. II are verified and hence the application 

zo ^ y^ 

which associates to (A~, Y_^) the fixed point of (|3.34p is continuously differentiable (in the sense of Frechet). 
When both A~ = and Y^ = 6 the Frechet derivative is the functional that maps {h^ , h^) G Zq into the 
function eA "^/i". 
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3.3.3 Analysis of the component of perturbation 

This step is devoted to the definition the component U^{t). First, we apply the change of variables introduced 
in Lemma [XH and we get that the matrix A{X^ ^ X^) in p.ip satisfies 

i(0, 0) = A" + Ar°, 

where and enjoy p.lPp and p.20p respectively. Relying on Remark [3131 we can still exploit estimate 

dSISD. 

We impose that X(t) = ^"(t) +y"*(r) + [/^'(r) is a solution of We then write C/P(r) = (C/^, U^f 

and, subtracting (|3.26p and p.33p from (|3.9p . we get 



A-{X- + t/-, + [/") - A-{X', 



u- 



X' 



dU^/dT = A"U" + N^U" + A°{x- + U-, X° + U°) - i*'(0, 0) 



(3.36) 



X'' 



Here, A" = ^-(0,0). 

Let be the metric space p.32p and consider the application Tp, defined for C/°) G 3^^^ as follows: 



T^i(C/^ f/°)[r] = 



r2(t/-, f/0)[r] 



+ [A-{X-{s) + U-{s), X°(.s) + C/°(s)) -A-(0, 0)]C/-(s)}ds 

e^"(--^){[7V"+i°(X-(s) + ;7-(s), X°(.s) + C/"(s)) -i"(0, 0)]c/°(s) 
i"(x-(s) + L/-(s), X"(s) + [/"(s)) -i°(0, X°(s))]x"(s)}ds 



(3.37) 



+ 



In the previous expression, X is the solution of (I3.33P and X^ is the solution of (|3.26p . We want to show 
that Tp maps into itself, provided that 6 is sufficiently small. We have 



C/°)[t]| < C^e-^'^^-^^[L\^U-{s)\ + \U\s)\ + \X\s)-Y^\]\X-{s)\ 
(s)| + \U-{s)\ + \X\s)\ + \U\s)\^ \U-{s)\ys 



< C. 



< 



k- iX-le-"'/^ + 2fcp<52e-(^+'^)/4 + 2S] kp5'^e-<''+''^l^ds 



c — a 
4 



C-Lkpk^S 



-LC-kp6^ 



^2g-r(3c+a)/4 



c + 2e 



LC_fcofc-(52e-^(2c-4s)/4^ 



" L C k kp S 



^2g-r(3c+a)/4 



8 



-6 



3c — a 



(3.38) 



In the previous expression, C_ is the same constant as in (13. 8p and L is a Lipschitz constant of yl~(X~, X") 
with respect to both the variables X~ and To obtain (|3.38p we exploit p.l5p . p.l6p . p.35p and the 
fact that, belonging to yg^, {U~ , C/°) satisfies 



|L/-(t)|, |[/0(r)| < fcp^2g-r(c+a)/4_ 

Also, the term ( is the same as in p.25p and we rely on the fact that |A~|, |C| < S. 



(3.39) 
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In the following expression L denotes a Lipschitz constant of A^{X , X^) with respect to both the 
variables X' and Also, we exploit the estimates ((XT^ . (lOH) . (jX^ . (jX^ and ([TT^ . 

|r^2(;7-, C/")[r]| < r |7VOC/0(s)|+i[|X-(s)| + |C/-(s)| + |XO(,s)| + |;70(s)|l|C/0(s)|ds 



4-00 



\X-{s)\ + \U-{s)\ + \U\s)\ \X\s)\ds 



+ 00 



< / — |f/°(s)|ds + L / [A:_|X-|e-^"/2 + 2fcp<52e-^(^+'^)/* + 2(5lfcp(52g-.(c+a)/4^g 



+ 00 



< 



/T /"T 
|X-|e~""/22Ms + L / 2fcp52e-^("+'^)/42<5ds 
,-00 ■/ +00 

Af(c+a) 3c + a c + a c + a 



-ct/2 



IQL5 



kp6 e 



-T{c+a)/A 



(3.40) 



Combining (I3.38|) and ()3.40p we get the following. Assume that the constant kp in p.32p is sufficiently large 
(namely, kp > ALk^/c). Then for every a < c — 4e we can choose 5 and M in such a way that Tp take values 
into y^^- Also, estimates similar to p.38p and (|3.40p ensure that one can choose the constants in such a way 
that Tp is a strict contraction. As a remark, we point out that, the bigger is a, the smaller is 5. 

We set a = 12e and we choose 6 in such a way that is a contraction from y^i2e itself. The constant 
e is the same as in p.l2p . However, in the following we regard as a map J^^q, where is the 

space p.32p obtained setting a = 0. In this way, we obtain that is a contraction on y^^, but, thanks to 
our choice of S, the fixed point automatically satisfies the sharper estimate 



|i7-(r)|, |C/°(r)| < kpS^e-^^^+^^''>/\ 
Also, in the definition of the space y^^ one can take = \(\ \X^ \ and hence 



(3.41) 



(3.42) 



where 2£. is defined by (j3.33p . Also, to simplify the notations in the previous expression we denote by C the 
point obtained applying the change of coordinates introduced in Lemma 13.21 to the vector (C, 0) defined by 



3.3.4 Prechet differentiability of the component of perturbation 

We are now concerned with the Frechet differentiability of the fixed point of the map Tp defined by (I3.37p . 
Since Tp{U~ , U^) depends on X~ and X'^, we regard Tp as a map 

T:yP,xy^yP,. (3.43) 

In the previous expression, y = y^ x y'~', where y~ and 3^" are defined by (|3.27p and p.28p respectively. 
Also, they satisfy £ 3^- and X° e y°. 

The proof of the differentiability relies on Lemma [3TT] (taking Y — y and X — y^f^)- We thus verify that 
the hypotheses of Lemma [01 are satisfied. 
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To simplify the exposition, we write (j3.37p as 



+ g(x-(.s) + U-{s), X"{s) + [/0(s))[/-(s)} 



ds 



+ l(x-(,s) + ;7-Gs), U\s))x^{s)]ds, 
where the functions F, G, H and L satisfy 

f(^x-{s), 0, x°{s), = 0(^0, o) EE h{{), = l(o, x", 



(3.44) 



0. 



(3.45) 



Note that X^{s) = r° is an equiUbrium for ^iTM) . 

Relying on (I3.41|) . one can show that the condition ()3.4|) is verified here, so applying Remark l3.1l we get 
that the fixed point (C/", U^) is Lipschitz continuous with respect to (X", X~). 

Concerning the Frechet differentiability of Tp with respect to [X^ , we proceed as follows. Fix 

an element (C/°, U' , X°, X") £ J^f^ x satisfying the estimates (IXT^ . (IX^ and (imjl . The 

Frechet differential of Tp with respect to X°) computed at the point ([/°, C/^, X^) is a linear 

map T e C{y, 3^0). The image of the element {h^ , hP) E y = x y" is given by 

Tp\h~, h°)[r]= J\^-(^-^^{f[x-{s), U-{s), X'^is), U"{s))h-{s) 
+ [Dx-F(^X-{s), U-{s), X°{s), U°{s))h-{s)\x-{s) 
+ [Dx-G(^X-{s) + [/-(,s), X°(,s) + [/°(,s))/i-(,s)] U-is) 
DxoF(^X-{s), U-is),X^{s), U'{s))h^{s)]x-{s) 
DxoG(^X-{s) + U-'{s), X°{s) + U°{s)'^h°{s) C/-(s)}ds 



Tp\h-, /i°)[r] 



^'^°'^^-'^{[Dx-H(^X-{s) + U-{s), X°(s) + C/"(s))/i-(s)] (/"(s) 



(3.46) 



+ [Dx-L[X-{.s) + U-{s), X^is), U%s)jh-{.s)\X%s) 
+ [DxoH(^X-{s) + U-{s), X°(s) + f/°(s))/i°(s)] U°{s) 
+ l(^X-{s) + U-{s), X°(s), U°{s)y°{s) 

Dx«l(x-{s) + U-{s), X"(s), t/"(s))/i°(s)]x°(s)}ds 



In the previous expression, [Dx-F{X-{s), U-{s),X°{s), U°{s))h-{s)] denotes the differential of the 
matrix valued function F with respect to the variable X^. The differential is computed at the point 
{X-{s), U-{s), X°{s), U°{s)) and is applied to the vector h-{s). To prove that indeed 



{T\h-,h°),T\h-,h^)) eyi 



one exploits estimate p.4ip and the identity L{0, 0) = 0. 

We now discuss the the Frechet differentiability of Tp with respect to ([/*', C/~). Fix an element 
{U°, X°, X-) e y'g^ X y. The Frechet differential of Tp with respect to (f/°, U'), evaluated at the 
point {U^, , X^, X^), is a linear map S G C{yQ, J^q) and the image of the element {h~ , /i°) G is 
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given by 



P-4 {t-s 



(3.47) 



+ g(x-(s) + U-{s), X"{s) + U"{s)y-{s) 

Du-G(x-{s) + U-{s), X'>is) + C/"(s))/i-(s)] U-is) 
DuoF(^X-{s), U-{s), X"{s), U'{s))h^{s)]x-{s) 
DuoG(^X-{s), U-{s), X°{s) + U"{s)'^h"{s) U-{s)^ds 
S^{U-, U°)[t] = e'^°^^-'^{[Du-H(^X-{s) + U-{s), X" (s) + U" {s)'^ {s)]u" {s) 

+ [Du-L(^X-{s) + U-{s), X°{s), C/"(s))/i-(s)]x"(s) 

+ N%°{s)+h(^X-{s), U-{s), X°{s) + U''{s)'^h°{s) 

Du«h(x-{s) + U-{s), X°{s) + U°is) 

DuoL[x-{s) + U-{s), C/''(s))/i°(s)]x°(s)}ds 

One can verify that, if ([/", , X", X^) G x y, then indeed S{h^ , h^) G y^. Also, S is continuous as 
a map from x 3^ in C{yl, y^). 

This shows that the hypotheses of Lemma [01 are all verified. 

3.3.5 Conclusion 

Applying Lemma |3. 11 we get that the map 

3^ - (3-48) 
that associates to (X^, X^) the fixed point of (|3.36p is Frechet differentiable and that its differential when 
X^{t) = and X^{t) = is the functional that associates to (/i^, h^) G y the functions U^{t) = 0, 
U^{t) = 0. We then perform the linear change of variables which is the inverse of the change of variables 
introduced in Lemma 13.21 In this way, we go back to the original variables. To simplify the notations, we 
still denote by {u^{t), U^{t)^ the functions obtained applying the change of variables. 

To define the map that parameterizes the uniformly stable manifold we proceed as follows: the orbit 
X^(t) is fixed. For every X_ G M"- , there exists a unique solution of p.33p . Also, in Section [5.3.2l we showed 
that the map 

X^X-{t) (3.49) 

is continuously differentiable in the sense of Frechet. As a consequence, the map obtained composing p.49p 
and p.48p is Frechet differentiable. Note that such a map associates to X~ the functions {X~ , U~, [7°). 
The function (j) that parameterizes the uniformly stable manifold is then defined by setting 

cl^{X)^[x-{0), X°(0) + C/"(0)). 

Thanks to the previous considerations, is continuously differentiable and the manifold is tangent to the 
stable space {{X~, 0) : X~ G M~} at the origin. Also, estimate (|3.24p is a consequence of (|3.42p . 
This concludes the proof of Theorem 13.11 



3.4 Uniformly stable manifolds 

Let be a fixed center manifold for the equation 

dU 



F{U), 



(3.50) 
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which satisfies Hypotheses [2] and [3] introduced in Section [21 In Theorem 13.11 we consider a fixed orbit lying 
on V'' and we construct the uniformly stable manifold relative to that orbit. In this section we discuss what 
happens if, instead of having a single orbit, we have a whole invariant manifold. 

More precisely, let Sq be an invariant manifold for p.9p and assume that Sq is entirely contained in the 
center manifold {X~ — 0}. Also, denote by S° the tangent space to Sq at the origin. Choosing a sufficiently 
small constant in Hypothesis [H we can assume that Sq is parameterized by S°. By construction, is 
contained in {X^ — 0}. Also, as in Section [3.1 .31 assume that Zq ~ {{X^ , 0, uq) : ( — 0} is a. manifold of 
zeroes for the function /o in (|3.6p . 

As a consequence of Theorem 13.11 we get the following result: 

Proposition 3.1. Let Hypotheses\^ and\^hold. Let Sq be an invariant manifold for p.9p entirely contained 
in the center manifold {X^ =0}. // the constant S in Hypothesis\^ is sufficiently small, then the following 
holds. 

There exists a continuously differentiable manifold which is defined in the ball of radius S and center 

at the origin. Also, satisfies the following properties: 

L A^^j^ is locally invariant for p.9p . meaning that if the initial datum lies on the manifold, then the 
solution (^{X^{t), Ar'^(r)) of p.9p lies on A^^^ for \t\ sufficiently small. 

2. A^5q is parameterized by S*' x and it is tangent to this space at the origin. Here, S*^ is the tangent 
space to at the origin and = {{X^ , 0) : X'^ = 0}. 

3. Any orbit Y{t) lying on A^^^ can be decomposed as 

Y{r) = r°(r) + Y- (r) + YP{t), (3.51) 

where Y'^{t) = (O, C^(t), uo{t)) is an orbit lying on Sq. The component Y^{t) ~ (X^(t), 0, O) lies 
on the stable manifold and the perturbation term YP{t) satisfies 

\YPiT)\ < C|C"(0)| |r-(0)|e-"^/^ (3.52) 

for some positive constant C . In p.52p . the constant c > is the same as in p.Sp . 

In the following we call A^^^ the uniformly stable manifold relative to Sq. 

Proof. Let (O, X"(t)) and (X~(r), O) be two orbits of (j3.9p lying on the center manifold {X^ — 0} and on 
the stable manifold respectively. We then have X'^{t) e y^, X~{t) G , where the metric spaces y^ and 
y^ are defined by (|3.30p . As in Section [X^ we use the notation y = y^ x 3^°. Consider the map 

'p-.y^yxyp, 

which associates to X-{t) and X°(t) the function (X-(r), X°(t), U-{t), U°{t)), where ([/", U°) is the 
perturbation term constructed in Section [3.3.3l We recall that J^fg is the set obtained setting a = in (I3.32p . 
As shown in Section fS. 3. 4[ the map $ is continuously differentiable in the sense of Frechet. Also, let 

/- : {^0 = 0} X {X- - 0} ^ y^ 

be the map that associates to (X^ , C, vP) G Zq the unique solution of the Cauchy problem (|3.33p . We recall 
that in p.33p Fq denotes (0, Uq). As shown in Section [3.3.21 the map is continuously differentiable in 
the sense of Frechet. Also, let 

f : {X- = 0} 3^,0 

be the map that associates to (X°, 0) the unique solution of the Cauchy problem ()3.26p . The map is also 
continuously differentiable in the sense of Frechet. Finally, fix a continuously differentiable map 

5° : S° ^ y° 
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parameterizing Sq. Define the map 
setting 

^{X", X-) = $ 1^/- (X-, g{X^)) , f o g%X"?j ■ (3.53) 

The map tp is then continuously differentiable in the sense of Frechet. By construction, V^fX". 2L^) is an 
element in the form (X~(r), X^{t), U^{t), U°{t)) and, setting 



Y{r) = [x-{r) + U-{r), X"(r) + [/"(r)) , 



we get that Y{t) can be decomposed as in (|3.51l) . Also, the perturbation term (U'~', U ) automatically 
satisfies p.52p . We then define the map 

parameterizing A^g^ by setting 

MX?, X-) = (x-{Q), x\Q) + u\Q)) = r(o), 

where X°(r), X-{t) and U°{t) are given by (^351) . 

The map -00 is continuously differentiable, being the composition of maps that are continuously differen- 
tiable in the sense of Frechet. Also, by construction the manifold TW^^ is invariant for p.9p . To prove that 
the manifold M.'g'^ is tangent to S" x V~ at the origin it is enough to observe that the Frechet differential 

of at = is the functional h_ ^ ^bT , while the Frechet differential of at X? = is the 
hmctional i-^ e"^°'^li' . □ 



4 Invariant manifolds for a singular ODE 

In Section [3] we extend to the general case the considerations introduced in Section 11.11 in the case of a toy 
model. In doing this, we apply the results obtained in Section [3] to study the singular ordinary differential 
equation 

Actually, most of the time we focus on system 

^=F{U). (4.2) 

We discuss several situations where (|4.ip and (|4.2p are equivalent, namely the Cauchy problem 

dr/dt = l/C[U{t)] 
r(0) = 

defines a continuously differentiable diffeomorphism r : [0, +oo[^ [0, +oo[. 

In Section [4] we exploit all Hypotheses [1] . . . [H Also, we rely on Proposition 14.11 whose proof is given in 
Section Before stating it, we have to introduce some notations. Let N denote the dimension of U. Also, 
?i_ is the number of eigenvalues of DF{0) with strictly negative real part, while (no + 1) is the number of 
eigenvalues of DF{0) with zero real part. Each eigenvalue is counted according to its multiplicity. Thanks 
to Hypothesis O N = + uq + 1. 
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Proposition 4.1. Let Hypotheses{l\ ... hold. If the constant 6 in Hypothesis\^ is sufficiently small, then 
in the ball with radius 6 and center at the origin we can define a continuously differentiable difjeomorphism 
T satisfying the following properties. Write T{U) — U as a column vector: 




where ( eR, uq e and m_ G M"- . // U satisfies (142]) then U satisfies 
dC/dr = Gio(C, Mo)uoC^ + Gi_(C, uq, u_)u_C 

< rfuo/rfT = {Goi(C, ^^o)+ [Go-(C, "0, ^^-)-Go-(C, uo, 0)]}c^io (4.3) 

du^/dr = Gs(C, Uq, u^)u^ 

In the previous expression, Gio is a row vector belonging to M"", Gi_ is a row vector in R"" , the matrices 
Goi and Go- belong to M""^"" and the matrix Gs belongs to M"-^"-. 

A center manifold of system ()4.3|) is the subspace {{(^, uq, 0) : U- = 0}, the stable manifold is the 
subspace {(0, 0, u_) : C — — 0}. Let M."^ be the uniformly stable manifold relative to the manifold 

E = {(C, 0, 0) : Uq = 0, U- — 0}, which is entirely constituted by equilibria. Then M'^ — {(C, 0, w_) : uq = 0}. 

In the statement of Proposition 14.11 by uniformly stable manifold relative to E we mean the manifold 
defined by Proposition l3.1l Also, note that by construction all the eigenvalues of the matrix Gs(0, 0, 0) have 
strictly negative real part. 



4.1 Slow and fast dynamics 

Let E denote, as before, the manifold of equilibria {(C, 0, 0) : uq = 0, u_ = 0}. 

Definition 4.1. A manifold of slow dynamics is a center manifold of (I4.3|l . In the following we fix the 
manifold of the slow dynamics {u_ = 0} and we denote it by A^". 

The manifold of fast dynamics of system (|4.3p is the uniformly stable manifold relative to E^ namely the 
subspace {uq = 0}. 

Note that both these manifolds are invariant for system (|4.3|) . Also, for every point (C, 0, u_) belonging 
to the manifold of fast dynamics, denote by (C('''), 0, u_(r)) the solution of (|4.3p such that 

(C(0), 0, u_(0)) = (C, 0, u-)- 

Combining (j3.5ip and (j3.52p we get that this solution decays exponentially fast to an equilibrium point. 
Namely, there exists Coo such that 

hm e'^'^\u-{T)\=Q= hm e^^/^KM - Cool, 

r — > + oo r — *+oo 

where the positive constant c satisfies Re\ < — c for every A eigenvalue of Gs(0, 0, 0). 
Consider system (j4.3p reduced on the manifold of slow dynamics: 

dC/dr = C^Gio(C, uo)uo 

duo / dr ^ GoiiC, uo)uoC (4.4) 

W_ EE 

If one goes back to the original variable t obtains 

dC/dt = CGio(C, uo)uo 

duo/dt = Goi(C, uo)uo (4.5) 

u- = 0, 
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namely an equation with no singularity. Note that (|4.4p and (14. 5p are equivalent. Indeed, by the uniqueness 
of the solution of a Cauchy problem associated to (|4.5p . the following holds. If C(0) > then ({t) > for 
every t. Thus, the Cauchy problem 

dr _ 1 
'dt ~ 



r(0) = 



(4.6) 



admits a global solution r : [0, +cxd[^ [0, +oo[ whose derivative is always different from 0. Thus, T{t) defines 
a change of variables and (|4.4p is equivalent to (|4.5p . 
One of our original goals is to study the solutions of 



dU_ 
Itt 



lying on a center manifold. Let Ad^'^ be a center manifold for 
Then A^^^ is a center manifold for 



around the equilibrium point (0, 0, 0). 



dc 

dt 
duQ 



= |Goi(C, Uq) + [Goi(C, Uo, U-) - Go-(C, Mo, 0)] juo 



(4.7) 



du^ 

IT 



-Gs(C, Mo, U-)U- 



We collect these results in the following 

Theorem 4.1. Assume that Hypotheses {1\ .. . [S| are satisfied. There exists an invariant center manifold 
j^oa system (14. 7p around the equilibrium point (0, (D, (D) which is contained in the manifold of the slow 
dynamics. In particular, equation ()4.7p restricted to A^^" is non singular and every solution satisfies the 
following property: if (^{0) > 0, then (^(t) > for every t. 

Remark 4.1. Hypothesis [8] ensures that the manifold {U : C{U) = 0} is invariant with respect to the slow 
dynamics. This hypothesis is not necessary to define an invariant center manifold Ad'^^ contained in the 
manifold of the slow dynamics. However, it is necessary if we want that (|4.4p is equivalent to (|4.5p . namely 
that the change of variables defined by (j4.6p is well defined. To see this, we can proceed as follows. 
Consider the equation 

Assume that one proceeds as in the proof of Lemma and exploits Hypotheses [T] . . . [7] but does not exploit 
Hypothesis [S] The system one eventually gets, restricted on the subspace {u_ =0}, is 



dC/dr = C5i(C, "0, 0) 
duo/dr = Goi(C, uo)uoC 

W_ EE 

where gi is the same function as in (|4.26p and satisfies 

5i(z, 0, 0) = Vz. 

Going back to the original variable t, (|4.8p becomes 

dC/dr = giiC, uq, 0) 
duo/dr = Goi(C, uo)uo 

M_ = 



(4.8) 



(4.9) 
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Thus, even if wc do not assume Hypothesis [51 the equation 

^-<-' 

restricted on the manifold of the slow dynamics {u_ = 0} is non singular. Also, one can define an invariant 
center manifold yVf'^^ which contains only slow dynamics. 

Note, however, that if Hypothesis [5] is not satisfied it may happen that for a solution U lying on Ai'^'^ 
C(C/(0)) > but C{U) touches in finite time. An example is the following. 

Consider the equation 

dui/dt = —U2 
du2/dt — 1*2(1 — U2) 
duz/dt — —u^/ui 

and set ^ 

((U)=ui F{U) = U1U2, uiul{l - U2), ~U3^ . 

Then Hypotheses [l] [3] . . . [7] are satisfied, but Hypothesis [8] is violated. The manifold of slow dynamics is the 
subspace {1*3 = 0} and it coincides with the center manifold A^"*'. Restrict to this subspace and consider 
the equation 

du2/dt — ^2(1 — U2). 

If < U2(0) < 1, then < it2(t) < 1 for every t. Also, du2/dt > for every t and hence ^2(0) < U2{t) < 1 
for every t. Since 

dui/dt = —U2, 

then by a comparison argument ui{t) < wi(0) — U2{0)t for every t > 0. In other words, if ui(0) > then 
ui{t) attains the value for some t < 'Ui(0)/u2(0). 

4.2 Applications of the uniformly stable manifold to the analysis of a singular 
ordinary differential equation 

We first recall a preliminary result we need in the following 

Lemma 4.1. Let C(r) be a real valued, continuous and bounded function satisfying ^(r) > for every 
T G [0, +cxd[. Let ^(t) be the maximal solution of the forward Cauchy problem 

dt/dr = Cir) .4 

Then t{T) is defined on the whole interval [0, 4-cxd[. Also, the following statements are equivalent: 
L t{T) is a continuously differentiable diffeomorphism t : [0, +cxd[^ [0, +cxd[. 

p + QO 

2. / C,{T)dT = +00. 



Condition 2 guarantees, in particular, that the inverse map T{t) is defined on the whole interval [0, +cxd[ 
and that it is continuously differentiable there. Also, note that C,{t) — (^{T{t)^ is automatically strictly bigger 
than for every t. 

Before stating the most important result in this section we need to introduce some notations. As before, 
c > denotes a positive constant satisfying ReX < —c for any A which is either an eigenvalue of Gs(0, 0, 0) 
or an eigenvalue with strictly negative real part of of Goi(0, 0). We denote by the subspace 

V'- = {(0, 0)}, 
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where ^ G M"" belongs to the eigenspace of Gio(0, 0) associated to the eigenvalues with strictly negative real 
part. Also, 

= {(0, C 0)}, 

where ^ e W^" belongs to the eigenspace of Gio(0, 0) associated to the eigenvalues with non positive real 
part. Clearly, 1/°^ C With V we denote the stable manifold: 

V— = {(0, 0, M_) : u- e M"-}, 

Finally, as in Section I^H we denote by E the manifold of equilibria {{(, 0, 0) : C G I^}- 
The most important result in this section is the following: 

Theorem 4.2. Let Hypotheses]^ ... [3 hold. If the constant 5 in Hypothesis]^ is sufficiently small, then in 
the ball with radius 5 and center at the origin one can define two manifolds, and tVC^", satisfying the 
following properties: 

1. both and Ai'^'^ are locally invariant for (14. 3p . namely: if the initial datum lies on the manifold, 
then the solution ((C(t), uo{t), w_(t)) of (|4.3p also lies on the manifold for \t\ sufficiently small. 

2. Ai" is contained in . 

3. AA^ is parameterized by E (B © V and it is tangent to this suhspace at the origin. Also, tVC^* is 
parameterized by E ® V^^~ © V and it is tangent to this subspace at the origin. 

4. let U{t) = (yC{'^)i uq{t), u_(t)^ be an orbit lying either on AA'^ or on AA'^^ and satisfying ^(0) > 0. 
Then the maximal solution of the forward Cauchy problem 

dt/dr = C(r) 

t{0)=0 ^ > 

defines a continuously differentiable diffeomorphism t : [0, +oo[^ [0, +oo[. Letrlt) denotes its inverse. 
Then the function U(t) ~ J7(T(i)) is a solution of ()4.H) and satisfies ({t) > for every t >0. 

5. any orbit lying on A4^ can be decomposed as 

U{t)^U-{t) + U'\t) + UP{t), (4.12) 

where U^{t) satisfies 

|C/-(t)| < fc_e-'=^/2|f/-(0)| (4.13) 

for a suitable constant fc_. Conversely, the component U^^{t) — (C''K^)j ''^oi''')^ ^ ^^^^ '"^ manifold 
of the slow dynamics. Also, if we use the variable t defined as the maximal solution of the Cauchy 
problem (j4.11l) . we have that the following property is satisfied. Denote by and Uq the first and the 
second component of U respectively. Then there exists a point (Coo, 0) such that 

hm (|C(t)-Coo| + |i^o(i)|)e^*/2=0. (4.14) 

t^-\-QO 

Finally, the perturbation term is small in the sense that 

|C/^(r)| <fcp|C'(0)||C/-(0)|e-^^/4 (4.15) 

for a suitable constant kp > 0. 

6. any orbit U{t) lying on tVC^* can be decomposed as 

U{t):^U-{t) + U''\t) + UP{t), (4.16) 

where U~ and UP satisfy \U-{t)\ < k^e-^^/^lU- {0)\ and \UP{t)\ < fcp|C'*'(0)||L/-(0)|e-=^/'' re- 
spectively. Here fc_ and kp denote the same constants as in ()4.13p and ()4.15p . The component 
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U''\t) ~ ((^^'(r), u*'(t), 0) lies on the manifold of the slow dynamics. More precisely, the follow- 
ing holds. Consider the maximal solution of the Cauchy problem 

dt/dT = C\t) 
t{Q) = 

and set C,'\t) = C^-rit)) and u'\t) = u''^{T{t)). Then (C^t), -u*'(t)) is a solution lying on a center- 
stable manifold of 

dC/dt = CGio(C, uo)uo 
dua/dt = G'oi(C, uo)uo 

W_ EE 0, 

Note that, strictly speaking, in (|4.12p and in (|4.16p the component docs not he on the manifold of 
the fast dynamics. Indeed, as we will see in the proof, is a solution of p.33p and hence docs not lie 
on {(0, 0, However, loosely speaking it can be regarded as a fast dynamic because of its exponential 

decay. 

Proof. We first define 

Consider system ()4.3|) restricted on the manifold of the slow dynamics. Thanks to the analysis in Section 
14.11 the variables t and r are then equivalent. Using the variable t, we get 

dC/dt = CGio(C, uo)wo 

duo/dt = Goi(C, uo)wo (4-17) 
u_ = 0, 

The manifold E = {{C, 0, 0) : C £ K} is then entirely constituted by equilibria. Applying Proposition 13.11 
to system (|4.17p with Sq — E, we then obtain MJ^*, the uniformly stable manifold relative to E, which is 
parameterized by © V'^^ . Note that so far we have used only the variable t: M]^* is a uniformly stable 
manifold for (|4.17p with respect to the variable t and by construction it is included in {u- = 0}, a center 
manifold for (|4.3p with respect to the variable r. The manifold is then obtained exploiting the variable 
T and applying Proposition 13. II to system (|4.3p with So = M"^^ . Also, the set 

Zo = {(0, 1*0, W-) : uo e M"o, u_ e M"-}. 

satisfies p.lOp . Properties 1, 3 and estimates (|4.13p and (I4.15P in the statement of Theorem 14.21 are then 
automatically satisfied, so we are left to prove estimate ()4.14p and property 4. 
To show that estimate (|4.14p holds we apply Lemma |4T] Thanks to (|4.12p . 

C(r)=C'(r) + C^(r), 

where [/^'(r) — (('"'(r), Uq{t), 0) lies on the manifold of the slow dynamics and C is the first component 
of the perturbation term W. Let t be defined as the maximal solution of 

di/dr = C^{t) 

m = 0, 

Then there exits (Coo, 0) such that 

Jim (IC'(f) - Cool + |ug'(^)l)e^*/' = 0. (4.18) 

Since |C^(t)| < kpd'^e^^'^^^, then for every r 

\iir)-tiT)\<Oil)S' 

where ^(t) is defined by (|4.1ip . Since also l^o (''') ^ "o(''')l — kpS^e^'^^^^, we conclude that (|4.14p implies 

Km . 
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Concering the proof of property 4, we apply Lemma |4. II Since W^ir) = {C\t), u'J{t), 0) lies on the 
manifold of the slow dynamics, then by the analysis in Section l4.1l it satisfies condition 1 in the statement 
of Lemma 14.11 and hence 



r + oo 

/ ('^ {T)dT = +00. 

Jo 



Since \Cp{t)\ < S^e-""^/^, then 

' iC' +e){r)dT^+oo. 
^ 

Applying again Lemma |4. II we get property 4. 

To define the manifold A^'^^ we proceed as follows. Consider M'^^, a center-stable manifold for (I4.17p . 

This manifold is parameterized by © V'^°~ and it is tangent to this space at the origin. The manifold Ai'^^ 

is defined applying Proposition 13.11 to system (|4.3p with Sq = Af^" and exploiting the presence of the set 

Zq = {(0, Mo, u_) : Mo e K"", u_ G M""} satisfying (|3.10p . Proceeding as before one gets that properties 

1, 3, 4 and 6 are satisfied. 

To verify property 2, we first observe that M"^" C M'^". To obtain Al" and A^'^" we applied Proposition l3.1l 
to iSo = M'^" and Sq = M'^^ respectively. Going back to the proof of Proposition 13.11 one can notice that 
the way we constructed the uniformly stable manifold with respect to is we associated to any orbit lying 
on Sq the manifold constructed in Theorem 13. II Thus the inclusion A/]^* C Af^* has as a consequence the 
inclusion C Al" . □ 



4.3 Proof of Proposition 14.11 
4.3.1 A preliminary result 

Before proving Proposition [5711 have to introduce a preliminary result, Lemma [ 

Let T be a continuously differentiable local diffeomorphism. To simplify the exposition, we also assume 
that T(0) = 0. Let := T{U) and 

F{U):^DT{T-^{IJ))F{T-^{U)) (4.19) 
If the function U{t) satisfies ()4.2p . then [/(r) solves 

^ = F{U)- (4.20) 

Also, given a real valued function Q{tj), let 

~aU):=Q[T-\U)\. (4.21) 
By direct check, one can verify that the following holds true. 

Lemma 4.2. Assume that Hypotheses\^\^ . . . El are satisfied by F and Also, assume that Hypothesis\^ is 

satisfied for some S. Then Hypotheses]^ ... are verified by F and C, and there exists 5, possibly smaller 
than S, such that Hypothesis]^ is as well satisfied. 



4.3.2 Proof of Proposition 14. It first part 

We are now ready to prove Proposition [13] The proof actually relies on standard techniques, but we give it 
for completeness. We proceed in several steps. 

• Step 1: let U = {ui . . . unY" be the components of U . Thanks to Hypothesis [H VC(0) ^ (J. Just to fix 
the ideas, we can assume 




By a smooth local change of variables we can assume that C(^) = ""i- Thanks to Lemma [4.21 Hy- 
potheses [T] . . . [S] are satisfied by the ODE written using the new variable. To simplify the exposition, 
we write U and C instead of U and C- 
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• Step 2: thanks to Hypothesis [51 there exists a manifold M^'' which is entirely constituted by equilibria 
and which is transversal to the manifold S, namely to {ui — 0}. Via a smooth local change of variables 
we can assume that the one-dimcnsional subspace 

E:={u2 = --- = UN = 0} (4.22) 

is entirely contained in A^*^'. Hypotheses [T] . . . [5] are satisfied in the new variables thanks to Lemma 
14:21 



• Step 3: let E be as in (|4.22p and denote by the eigenspace of DF{6) associated to eigenvalues 
with real part. Also, let V be the eigenspace associated to eigenvalues with strictly negative 
real part. The dimension of V and of V~~ is no + 1 and n_ respectively. Thanks to Hypothesis [31 
A'' = no + 1 + n_. The vector (1, ... 0) belongs to because E C V^. Also, we can assume, via a 
linear change of variables, that 

V = {Uno+2 = ...UN ^0} = {C = 0, U2 = . . . U„o + l = 0}. 

Fix a center manifold Ai'^ for system 

^ = FiU) (4.23) 

around the equilibrium point 0: M'^ is parameterized by and it is tangent to this space at the 
origin 0. Also, let M^jf be the uniformly stable manifold of system ()4.23p relative to the manifold of 
equilibria E defined by ()4.22|) : this manifold is paramerized hy V ® E and it is tangent to this space 
at the origin. By a local smooth change of variables we can assume that actually 

M" = {Un„+2 = ...UN = 0} Me' = {u2 = • ■ • Uno + l = 0}. 

Note that the Hypotheses [T] . . . [H] are satisfied because of Lemma 14.21 

• Step 4- consider the following decomposition: 

/ C \ / /i(C,"o, \ 

U=\ uo F{U)= Fo{C,uo,u^) (4.24) 

V / V F^iC, uo, u_), J 

where j\ € M, uq, Fq G M"" and m_, F_ G M."- . In the new coordinates, the center manifold A^'^ 
is the subspace {u_ = 0} and the uniformly stable manifold A4'^ is {uq = 0}. 

The center manifold {u_ = 0}is invariant for the equation 

^ = F{U) (4.25) 

and hence F-(C, uq, 0) = for every C and uq. By regularity, 

F-iC, Uo, u_) GsiC, Uo, u^)u^ 

for a suitable matrix Gg £ M"-^"-. Also, the uniformly stable manifold is invariant and hence 
proceeding as before we get that 

Fo(C, Uo, -"-) = Gc(C, Uo, u^)uo 
for a suitable matrix Gc G M"'^^'^"". Finally, Hypothesis [7] imphes that 

/i(0, Uo, U-) = 

and hence by regularity /i(C, uo, uJ) = 5i(C, uo, u^)Q. Consider the decomposition 
Ge(C, uo, u_) = Ge(C, uo, 0) + [Ge(C, uo, u_) - Gc(C, uo, 0)] . 
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Thanks to Hypothesis [5l the subspace {( — 0, — 0} is entirely constituted by equihbria and hence 

Ge(0, uo, 0) = 0. 

By regularity, GdC, uq, 0) = G'oi(C, mo)C for a suitable matrix Gqi G M""^"" . Putting all the previous 
considerations together, we get that system ()4.25p can be written as 



duo/dT = {Goi(C, uo)C + [Ge(C, uo, u^) - G,(C, uo, 0)] } 
du^/dr = Gs{Ci "Oj u^)u- 



Uo 



(4.26) 



Consider the decomposition 

5i(C, Mo, U-) = gi{C, Uo, 0) + [gi(C, "o, u-) - .gi(C, "o, 0)] 

By construction Gs(0, 0, 0) admits only eigenvalues with strictly negative real part, thus Gs{C,, uo, u^)u^ 
implies m_ = 0. Thus, the set {U : C{U) = 0, F{U) = 0} is the subspace {( = 0, u_ = 0}. Thanks to 
Hypothesis m we have 

31 (0, uo, 0) =0. 

By regularity, we thus have 

5i(C, Uo, 0) = 5ii(C, uo)C [giiC, Uo, u_) - gi{C, uo, 0)] = Gi_(C, uo, U-)u- 

for a suitable row vector Gi_(C, mq, m_) £ W"- . Also, since the manifold {uo = 0, u- = 0} is entirely 
constituted by equilibria, then gii{C, 0) = for every C and hence 

ffii(C, "o) = Gio(C, uo)uo 
for a suitable vector Gio £ R"". In other words, (I4.26P reduces to 



dC/dr = C^Gio(C, uo)uo + CGi-(C, uo, U-)u^ 

duo/dr = |Goi(C, uo)C + [Gc(C, "o, U-) - GdC, uo, 0)] juo 

du^/dr — Gs(C, Wo, U-)u- 



(4.27) 



;Step 5.' we introduce a refined change of variables. Consider system Ij4.26p restricted on the invariant 
subspace {C = 0}. One obtains 



duo/dr = Gc(0, uo, u_) - Gc(0, uo, 0) 
du-fdr — Gs(0, uq, U-)u- 



uo 



(4.28) 



The subspace {u_ = 0} is entirely constituted by equihbria. Also, given a point {uo, U-) belonging to a 
small enough neighbourhood of 0, then the solution of (|4.28p starting at {uq, u_) decays exponentially 
fast to a point in the subspace {u- — 0}. This is a consequence of the fact that Gs(0, 0, 0) admits 
only eigenvalues with strictly negative real part. 

We can define a change of variables U = T^([/) such that in the new variables U the following holds. For 
every uo{0) S M"" and for every ?I_(0) £ M"- , the solution of starting at the point (uo(0), U-{0)) 

converges exponentially fast to the point (uo(0), 6). In other words, the set {uq = mo(0)} is the stable 



33 



manifold of system (|4.28p around the equilibrimii point (uo(0), O)- Let F{U) be defined as in (|4.19p . 
with T = T*. Then 

(C^GiolC uo)uq + CGi_(C, uq, u_)u^ 
{Goi(C, ua)C+ [Ge(C "0, u^) ~ Gc(C, uo, 0)] }so 
Gs(C, Wo, -u_)?2_ 

Because of the way we chose T4, when C = then duo/dT = and hence 

[Ge(0, Uo, U-) - Ge(0, So, 0)]f2o = 0. 

By regularity, 

[Ge(C, UO, U-) - GciC 0)] = [Go-(C, "0, - Go-(C, ^^o, 0)]C 



for a suitable function Go- G M"" 



xno 



• Step 6: to conclude the proof of Proposition 14.11 it is enough to define the local diffeomorphism T as 
the composition of all the local diffcomorphisms defined in the previous steps. 
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